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In the stochastic mean-field approach, an ensemble of initial conditions is considered to incorporate 
correlations beyond the mean-field. Then each starting pont is propagated separately using the 
Time-Dependent Hartree-Fock equation of motion. This approach provides a rather simple tool to 
better describe fluctuations compared to the standard TDHF. Several illustrations are presented 
showing that this theory can be rather effective to treat the dynamics close to a quantum phase 
transition. Applications to fusion and transfer reactions demonstrate the great improvement in the 
description of mass dispersion. 



I. INTRODUCTION 



II. THE STOCHASTIC MEAN-FIELD THEORY 
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The mean-field description of a many-body system, i.e. 
the Hartree-Fock (HF) and/or time-dependent Hartree- 
Fock theory (TDHF), provides a simple tool for descrip- 
tions of certain aspects of complex quantum systems. 
However, it is well known that the mean- field approxi- 
mation is suitable for the description of mean values of 
one-body observables, while quantum fluctuations of col- 
lective variables are severely underestimated. A second 
limitation of mean-field dynamics is that it cannot de- 
scribe spontaneous symmetry breaking during dynamical 
evolution. If certain symmetries are present in the initial 
state, these symmetries are preserved during the evolu- 
tion. We have recently shown that a stochastic mean- 
field (SMF) approach fl|, @] where the TDHF evolution 
is replaced by a set of mean-field evolution with prop- 
erly chosen initial conditions. It will be shown that this 
approach can be a suitable tool to go beyond mean-field 
and describe the evolution of a system close to a quan- 
tum phase-transition Q. In a series of article, we ap- 
plied the SMF approach to describe transport properties 
in fusion reaction. Transport coefficients related to dis- 
sipation and fluctuations have been obtained that 
are crucial to understand the physics of Heavy-Ion colli- 
sions around the Coulomb barrier. A summary of recent 
results is presented. 



In a mean-field approach, the nuclear many-body dy- 
namical problem is replaced by a system of particles in- 
teracting through a common self-consistent mean-field. 
Then, the information on the system is contained in the 
one-body density matrix p that evolves according to the 
so-called TDHF equation: 
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where h[p] = d£(p)/dp denotes the mean- field Hamilto- 
nian. While quite successful in the description of some 
aspects of nuclear structure and reactions [7| , it is known 
to not properly describe fluctuations of one-body degrees 
of freedom, i.e. correlations. Numerous approaches have 
been proposed either deterministic or stochastic to ex- 
tended mean-field and describe fluctuations in collective 
space (see ref. [l| and reference therein). Most often, 
these approaches are too complex to be applied in re- 
alistic situations with actual computational power. A 
second limitation of mean-field dynamics is that it can 
not describe spontaneous symmetry breaking during dy- 
namical evolution. If certain symmetries are present in 
the initial state, these symmetries are preserved during 
the evolution Hi]. 

The Stochastic Mean-Field (SMF) has been recently 
shown to provide a suitable answer for the description of 
fluctuations as well as of the symmetry breaking process 
while keeping the attractive aspects of mean-field. Let 
us assume that the aim is to improve the description of 
a system that, at the mean- field level and time t , is 
described by a density of the form: 
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Note that, this density can describe either a pure Slater 
determinant (rij = 0, 1) or more generally an initial 
many-body density of the form: 

D = - exp ^i a \ a i) ( 3 ) 

where z is a normalization factor while (a|,aj) are the 
creation/annihilation operators associated to the canon- 
ical basis \ fi). Then, the mean- field evolution, Eq. ([T]) 
reduces to the evolution of the set of single-particle states 

d 

ih—\cpi(t)} = h\p]\tpi(t)), (4) 

while keeping the occupation numbers constant. 

In the SMF approach, a set of initial one-body densities 

p x (t Q ) = 5^|Vi>/&(«o)tol (5) 

ij 

is considered, where A denotes a given initial density. 
The density matrix components p* ■ are chosen in such a 
way that initially, the density obtained by averaging over 
different initial conditions identifies the density ©. 

It was shown in ref. that a convenient choice for 
the statistical properties of the initial sampling is 

Pi, j (to) = S ij n i + Spf j (t ), (6) 

where #py(to) are mean-zero random Gaussian numbers 
while 

H(*o)^r(*o) = I s * 6 * (<(! - 4) + »i C 1 - <x)) 

The average is taken here on initial conditions. In this ap- 
proach, each initial condition given by Eq. (J5j> is evolved 
with its own mean-field independently from the other tra- 
jectories, i.e. 

ih^{t)) = h\p x ]\^{t)), (8) 

while keeping the density matrix components constant. 
Therefore, the evolution along each trajectory is similar 
to standard mean-field propagation and can be imple- 
ment with existing codes. A schematic illustration of the 
standard mean-field and stochastic mean-field is given in 
figure [TJ 

III. FROM DENSITIES TO OBSERVABLES: 
EHRENFEST FORMULATION OF THE 
STOCHASTIC MEAN-FIELD THEORY 

The mean-field theory is a quantal approach and even 
if it usually underestimates fluctuations of collective ob- 
servables in the nuclear physics context, these fluctua- 
tions are non-zero. Within mean-field theory, the ex- 
pectation value of an observable O is obtained through 




FIG. 1: Illustration of the quantal mean-field approach (top) 
where a single density is propagated and of the SMF approach 
(bottom) where a set of densities are chosen and where each 
initial density is propagated independently from the others. 

(O) = Tr (OD) where D has the form ([3]). Accordingly, 
the quantal average and fluctuations of a one-body ob- 
servable Q along the mean-field trajectory are given by: 

<<9>=X>i(*)IQl¥>i(*)>ni (9) 
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and 

"h® = (Q 2 ) - (Q? = E \{vi(t)\Q\v>i(t))\ 2 m(i ~ ndio) 
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An important aspect of the SMF approach is that the 
quantum expectation value is replaced by a classical sta- 
tistical average over the initial conditions. Denoting by 
Q x (t) the value of the observable at time t for a given 
event, fluctuations are obtained using 

^mf(g = WW~oW (ii) 

where Q(t) — Q x (t). The statistical properties of ini- 
tial conditions insures that quantal fluctuations [oq] and 
statistical [og MF ] fluctuations are equal at initial time. 
Note that such a classical mapping is a known technique 
to simulate quantum objects and might even be exact in 
some cases [T(J EH ■ 

In practice, it might be advantageous to select few col- 
lective degrees of freedom instead of the full one-body 
density matrix. At the mean-field level, the evolution of 
a set of one-body observable Qi is given by the Ehrenfest 
theorem: 

ihj t (Q i ) = ([Q i ,H}) (12) 

If a complete set of one-body observables is taken, for 
instance if we consider full set of operators{a]aj}, one 
recovers eq. (fTJ). In many situations, one might further 
reduce the evolution to a restricted set of relevant degrees 
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of freedoms in such a way that the mean-field approxi- 
mation leads to a closed set of equations between them, 
i.e. 



ih-(Qi) = F [(Q x ) 



' i (Qn, 



(13) 



Starting from this equation, one can also formulate the 
SMF theory directly in the selected space of degrees 
of freedom by considering a set of initial conditions 
{Qi (io)}i=i,n an d by using directly the evolution: 



I Q?(i), ■•■,<&(*) 



(14) 



for each initial condition A. Note that statistical proper- 
ties, i.e. first and second moments, of initial conditions 
should be computed using the conditions (JSJ) and ([7]). 



IV. ILLUSTRATIONS 

In recent years, we have applied the SMF approach ei- 
ther to schematic models or to realistic situations encoun- 
tered in nuclear reactions where mean- field alone was un- 
able to provide a suitable answer. Some examples are 
briefly discussed below. 



A. Many-body dynamics near a saddle point 

As mentioned in the introduction, the mean-field the- 
ory alone cannot break a symmetry by itself. The sym- 
metry breaking can often be regarded as the presence 
of a saddle point in a collective space while the absence 
of symmetry breaking in mean-field just means that the 
system will stay at the top of the saddle if it is left here 
initially. Such situation is well illustrated in the Lipkin- 
Meshkov-Glick model. This model consists of N parti- 
cles distributed in two N-fold degenerated single-particle 
states separated by an energy e. The associated Hamil- 
tonian is given by (taking U = 1), 



H = eJ z - V(Jl - Jl 



(15) 



where V denotes the interaction strength while Ji (i = x, 
y, z), are the quasi-spin operators defined as 
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;(J + + J-), Jy = -{J+-J-) (16) 



with J + = Y^>=i c +.p c --p> J- = J+ an d where c+ p and 
c_ p are creation operators associated with the upper and 
lower single-particle levels. In the following, energies and 
times are given in e and h/e units respectively. 



It could be shown that the TDHF dynamic can be 
recast as a set of coupled equations between the expec- 
tation values of the quasi-spin operators ji = (Ji)/N (for 
i = x, y and z) given by: 
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where x — V(N — l)/e. Note that, this equation of mo- 
tion is nothing but a special case of eq. (fT3| where the 
information is contained in the three quasi-spin compo- 
nents. To illustrate the symmetry breaking in this model 
it is convenient to display the Hartree-Fock energy £hf 
as a function of the j z component (Fig. [5]). Note that, 
here the order parameter a = iacos(— j z /2) is used for 
conveniency. When the strength parameter is larger than 
a critical value (x > 1), the parity symmetry is broken 
in a direction. For x > 1, if the system is initially at the 
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FIG. 2: (color online). Evolution of the Hartree-Fock energy 
£hf as a function of a for x = 0-5 (dashed line), x — 1-8 
(doted line) and X = 5 (solid line) for iV = 40 particles. 
The arrow indicates the initial condition used in the SMF 
dynamics. 

position indicated by the arrow in Fig. [2 with TDHF it 
will remain at this point, i.e. this initial condition is a 
stationary solution of Eq. (IT71) . 

Following the strategy discussed above, a SMF ap- 
proach can be directly formulated in collective space 
where initial random conditions for the spin components 
are taken. Starting from the statistical properties (|6]) 
and [[7]). it could be shown that the quasi-spins should 
be initially sampled according to Gaussian probabilities 
with first moments given by |3(: 



j£(*o)=#(io) = 0, 
and second moments determined by, 



(18) 
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while the z component is a non fluctuating quantity. 



(19) 
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FIG. 3: (color online) Top: illustration of the initial sampling 
used for the SMF theory in the collective space of quasi-spins. 
Bottom: Exact evolution of the z quasi-spin component ob- 
tained when the initial state is \ j, —j) for three different values 
°f X- X — 0-5 (solid line), \ = 1-8 (dotted line) and x = 5-0 
(dashed line) for N — 40 particles. The corresponding results 
obtained with the SMF simulations are shown with circles, 
squares and triangles respectively, (adapted from Q]) 



An illustration of the initial sampling (top) and of re- 
sults obtained by averaging mean-field trajectories with 
different initial conditions is shown in Fig. [3] and com- 
pared to the exact dynamic. As we can see from the 
figure, while the original mean-field gives constant quasi- 
spins as a function of time, the SMF approach greatly im- 
proves the dynamics and follows the exact evolution up to 
a certain time that depends on the interaction strength. 
As shown in Fig. [4] the stochastic approach not only 
improves the description of the mean-value of one-body 
observables but also the fluctuations. 



V. APPLICATION TO NUCLEAR REACTIONS 

The SMF has been recently used to deduce transport 
coefficients associated to momentum dissipation or mass 
transfer during reactions from a fully microscopic the- 
ory 041] ■ The TDHF theory provides a powerful way 
to get insight nuclear reaction and treat various effects 
like deformation, nucleon transfer, fusion, ... in a quan- 
ta! transport theory. An illustration of nuclear densities 
obtained at various time of the 40 Ca + 90 Zr reactions is 
given in Fig. [5j 
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FIG. 4: (color online) Exact evolution of dispersions of quasi- 
spin operators obtained when the initial state is \j, —j) for 
three different values of x> from top to bottom \ = 0-5 (a), 
X — 1.8 (b) and x = 5.0 (c) are shown. In each case, solid, 
dashed and dotted lines correspond to cr x (t), <Jy(t) and A^(t), 
respectively. In each case, results of the SMF simulations 
are shown with triangles (<r^), squares (<jy) and circles (o~z)- 
(taken from [3J) 



The mean-field approach does include the so-called 
one-body dissipation associated to the deformation of 
the system and/or to the exchange of particles. For in- 
stance, considering a set of observables, denoted gener- 
ically Q = {Qi}, like the relative distance, relative mo- 
mentum, angular momentum between nuclei or the num- 
ber of nucleons inside one of the nucleus, it is possible 
to reduce the TDHF evolution and obtain classical equa- 
tions of motion of the form: 



dQ l 
dt 



F(Q,t)-5^i/«(Q,i)Qj 



(20) 



where F is an eventual driving force while v corresponds 
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FIG. 5: (color online). Snapshots of the nucleon density pro- 
files on the reaction plane, p(x,y,z = 0), are indicated by 
contour plots for the central collision of 40 Ca + 90 Zr system 
at i?cm = 97 MeV in units of fm -3 . The black dot is the cen- 
ter of mass point. The red lines indicate the positions of the 
window xo and Do = dxo /dt denotes velocity of the window, 
(taken from [(|) 



to drift coefficients. For instance, the nucleus- nucleus in- 
teraction potential and energy loss associated to internal 
dissipation has been extracted in ref. 0, [l2j using such 
formula. 

When TDHF is extended to incorporate initial fluctu- 
ations, the equation of motion itself becomes a stochastic 
process: 



d0[ 
dt 



F(Q\t)-J2^(Q\t)Qj+SQt (21 



For short time, the average drifts should identify with 
the TDHF one while the extra term is a random variables 
that leads to dispersion around the mean trajectory. In 
the Markov limit, one can define the diffusion coefficient 



28{t-t')D ij {t). 



(22) 



This mapping has been recently used to not only study 
dissipative process but also estimate fluctuations prop- 
erties in the momentum and mass exchange. Denoting 



by DAA{t) the diffusion coefficient associated with mass, 
fluctuations in mass of the target and/or projectile can 
be computed using the simple formula 



2 / D AA {s)ds. 



(23) 



In figure [51 an example of estimated variances during 
the asymmetric reaction 40 Ca + 90 Zr is shown as a func- 
tion of time in the case of fusion reaction (top) or below 
the Coulomb barrier (middle and bottom panel). All 



10 
8 
6 
4 
2 




< 

OJ < 





• ^exc 


• 


- SMF 


*/ 




*/ 






1 E cm =110MeV 






(a)" 
i , , , , 



t(10" 22 s) 



10 



< 

OJ < 



< 

OJ < 



10 
8 
6 
4 
2 


0.8 
0.6 
0.4 
0.2 




• ^exc 
- SMF 












- E cm =97 MeV 








■ iii 


(b)- 



10 



15 



20 



t(10" 22 s) 



1 1 1 1 1 1 1 r 

- • - N exc 


i _ 


' SMF 




m- 




- E cm =93MeV f 




# ~r"T I 1 I I i_ 


(0- 
. I . . 


5 10 


15 


t(10" 22 s) 





FIG. 6: Variances of fragment mass distributions are plotted 
versus time in collisions of 40 Ca + 90 Zr system at three dif- 
ferent center-of-mass energies. The dotted lines denote total 
number of exchanged nucleons until a given time t. (taken 
from [6(]) 

cases correspond to central collisions. Note that below 
the Coulomb barrier the target and projectile re-separate 
after having exchanged few nucleons corresponding to 
transfer reactions. In general, it is observed that the 
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fluctuations are greatly increased compared to the orig- 
inal TDHF and are compatible with the net number of 
exchanged nucleons from one nucleus to the other (dotted 
line in Fig. [5]). 

VI. SUMMARY 

In this contribution, illustrations of the application 
of the stochastic mean-field theory are discussed. It is 
shown, that the introduction of initial fluctuations fol- 



lowed by a set of independent mean-field trajectories 
greatly improves the original mean-field picture. In par- 
ticular, it seems that this approach is a powerful to in- 
crease the fluctuations that are generally strongly under- 
estimated in TDHF or to describe the many-body dy- 
namics close to a saddle point. 
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